The Fundamental Equations of Electron Motion 
(Dynamics of High Speed Particles) 

By L. A. MacColl 

I. Introduction 

In work relating to the motion of electrons and other particles it is fairly 
common to assume that the particles obey the laws of Newtonian dynamics. 
That is, briefly, it is assumed that the rectangular coordinates (x, y, z) of 
the particle under consideration satisfy the differential equations 

mx = Xj my = Y, mi = Z, 

where m is the mass of the particle (assumed constant), X, Y, and Z are the 
components of the applied force, and the dots indicate differentiation with 
respect to the time t. 

However, it is well recognized now that the above equations are not 
strictly correct, and that they merely represent an approximation which is 
adequate when the speed of the particle is sufficiently small compared with 
the speed of light. The system of dynamics based upon the correct equa- 
tions^ (which will be exhibited presently) is commonly called relalivisiic 
dynamics, not because any knowledge of the theory of relativity is essential 
to its understanding and use^ but because it is in agreement with the theory 
of relativity (which Newtonian dynamics is not), because it was first de- 
veloped in connection with work on the theory of relativity, and because 
even yet virtually all of the expositions of the subject are to be found in 
books and papers dealing primarily with the theory of relativity. 

Just where the dividing line should be set between cases in which New- 
tonian dynamics is an adequate approximation and cases in which it is 
necessary to use relativistic dynamics is, of course, a rather vague question 
which cannot be answered simply and definitely. We may note, however, 

1 It is not the purpose of this article to discuss questions of fundamental physics, or 
the physical validity of any particular equations. For purposes of discussion, we assume 
outright that relativistic dynamics is at least more nearly correct than is Newtonian 
dynamics. 

2 The theory of relativity can be described briefly as a theory of the relations between 
the descriptions of phenomena in terms of different systems of reference. We shall not 
be concerned with this theory, because we shall be employing the same reference system 
throughout most of our discussion. In the final section of the paper we shall consider 
purely geometrical transformations of the coordinate system. These transformations, 
however, involve nothing that is reaCy characteristic of the theory of relativity in the 
usual sense. 
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that according to relativistic dynamics the mass of a five thousand volt 
electron is about one per cent greater than the mass of an electron at rest. 
From this we can infer that, while Newtonian dynamics may be adequate 
for many purposes in our studies of electron motion, we do not have any 
great amount of margin, and that it wiU be necessary to use relativistic 
dynamics whenever we wish to obtam really good results concerning the 
motion of even moderately high speed electrons. 

This article is purely expositor}'. Its purpose is to set forth the funda- 
mental equations and theorems of relativistic particle dynamics in a clear 
and concise form, unencumbered with any material relating to the theory of 
relativity proper. Almost all of the material is to be regarded as already 
known, but apparently it is only to be found in an inconvenient and scattered 
form. The incomplete bibliography at the end of the paper gives references 
to some of the more accessible sources of this and other related material. 

II. The Elementary Differential Equations of Motion 

Our discussion might be begun in any one of a number of ways, and no 
doubt the different approaches would appeal unequally to different readers. 
Considering the nature and purposes of this article, the author has deemed 
it best to begin by writing down at once the differential equations of motion 
of a particle (according to relativistic dynamics) in their most elementary 
form. Then, for the purposes of this discussion, these equations will have 
the status of a fundamental assumption. It need hardly be said that the 
equations are not written down arbitrarily. On the contrary, they represent 
the consensus of modern opinion as to the laws under which particles really 
do move.^ The grounds, experimental and theoretical, for this opinion are 
set forth in various of the works cited in the bibliography. 

For the time being, until the contrary is stated in the final section, we 
employ a fixed rectangular coordinate system. Instead of denoting the 
coordinates of the particle by x, y, and s, as we have done provisionally in 
the Introduction, we shall denote them by .ri, .T2, and .T3. Then ±1, ±2, and 
X3 denote the components of the velocity of the particle. The components 
of the force acting on the particle will be denoted by .Yi, X2, and X3. For 
the time being we need only note that the force may depend upon the 
coordinates, the velocity, and the time; later on we shaU introduce some 
more explicit assumptions about the force. The symbol c will be used to 
denote the speed of light in vacuo. 

" The validitj' of these laws is not unrestricted. It is hmited on the one hand by the 
quantum phenomena which become appreciable on the atomic scale, and on the other hand 
by certain phenomena revealed by the general theory of relativity which become 
appreciable on the cosmic scale. 
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We assume that the particle moves, under the influence of the force (Xi, 
.Y2, Xi), so that its coordinates satisfy the system of differential equations 



ttlfiXn 



dt Vi - (z'-A-) 



^ X„, (■« = 1,2,3), (1) 



where m» is a positive constant characteristic of the particle, and v- is an 
abbreviation for the expression ±1' + xi' + ±3^* The positive value of the 
square root is the significant one; and wherever square roots appear in the 
subsequent work it will be understood, unless the contrary is stated, that 
the positive values are intended. 

A few remarks may help bring out the significance of the foregoing assump- 
tion and its relations to the corresponding fundamental assumption of 
Newtonian dynamics. 

We call the constant mo the rcsl-mass of the particle, and we assume {in 
accordance with the experimental evidence) that im is identical with the 
mass of the particle which is used in Newtonian dynamics. In relativistic 
dynamics the quantity m defined by the equation 



ma 

m = 



Vl - (wVc'} 

is called the mass of the particle. We note that as v/c approaches zero the 
mass approaches the rest-mass (whence the appropriateness of the latter 
term), and that as v/c approaches unity the mass increases without limit. 
Consider the vector having the components pi, pi, ps defined by the 
formulae 

We call this vector the momentum of the particle. The momentum is equal 
to the velocity of the particle multiplied by the mass. 

Now equations (1) assert that the time-rate of change of the momentum 
of the particle is equal to the applied force. 

We have already obser\ed that as v/c approaches zero the relativistic mass 
of a particle approaches the Newtonian mass. We now note that as v/c 
approaches zero the components of the relativistic momentum approach 
the values 

p„ ^ Wo in, (2') 

* We might merely say that v is the speed of the particle. However, for our immediate 
purposes, it is important not to lose sight of the fact that i- is a certain particular function 
of the components of velocity. 
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which are precisely the components of the momentum according to the 
Newtonian theory. 

Finally, as v/c approaches zero, the difierential equations of motion (1) 
approach the forms* 



^ (mox„) = Xn, 



(10 



which are the Newtonian differential equations of motion. 

Thus we see that Newtonian dynamics is in effect a simplified approximate 
form of relativistic dynamics which is valid when the speed of the particle 
under consideration is sufficiently small compared with the speed of light. 

Let us carry out the indicated differentiations in equations (I), and then 
solve the resulting equations for the quantities JUoXi, moXn, nioXs. The work 
is straightforward, and need not be given here. We obtain the following set 
of formulae: 



2 -2\-l/2 



/< 2 — 2\ 

WuXi = (1 — If c J 



2 _-2i-l/2 



WflXa = (1 — ^ C ) 



2 -2^-l/^ 



Xi 


Xi i-2 c 




XiXiC 




X2 


l-(il^ + X3= 


)c-'' 


±2 ±3 C^^ 


> 


X, 


±2 Xz C~^ 




\~{x,' + x,%-' 




1 - 


ixi + xi)c-' 


Xi 


Xi x^ c~^ 






Xi ±2 C~' 


X2 


X2 Xz c~^ 


, (3) . 




Xi ±3 C~^ 


Xz 


1 - (Xi' + X2')C-' 




1 - 


{Xi i- X3 )c 




Xi ±2 C"^ ^1 






Xi ±2 C~ 


1 - 


(Xi' + X3')c-' X2 


. 




±1 ±3 C~^ 




X2X3C~^ Xz 





ntoXz = {I — V c 



These equations are, of course, the differential equations of motion (1) 
written in a new, but equivalent, form. 

If, at some particular instant, the particle is moving parallel to the .Vi-axis, 
so that X2 = Xs ^ 0, tlie equations (3) reduce at that instatit to the forms: 



ntoXi 



= X,, 



misX2 



= X2, 



inoXz 



= Xz. 



These equations show that a particle of rest -mass mo, moving with speed v, 
responds to a force parallel to the velocity as would a Newtonian particle* 
of mass 

^l - ("l"- ^2^2)3/2'' 

* If this conclusion is not entirely evident, the reader is referred to equations (3), from 
which the conclusion follows at once. 

' I.e. an ideal particle which obeys the laws of Newtonian dynamics. 
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and that the particle responds to a force perpendicular to the velocity as 
would a Newtonian particle of mass 

mo 

mt = 



(1 - v'(r^Y''-' 

For this reason, it was usual in the early work on relativistic dynamics to 
ascribe two masses to a particle: the longitudinal mass m., and the transverse 
mass mt. However, in general this procedure leads only to inconveniences, 
and it has been almost entirely abandoned. 

This concludes our discussion of the elementary differential equations of 
motion. Without any further general theory of relativistic dynamics it is 
possible to solve many interesting and important problems. For instance, 
it can be shown easily that the trajectory of a particle subjected to a force 
which is constant in magnitude and direction is a catenary (rather than a 
parabola, which is the curve predicted by Newtonian dynamics).^ In the 
following sections we shall discuss some of the less elementary parts of the 
subject. 

III. The Lagrangian Equations 

In the foregoing the components of the applied force have been any func- 
tions of the coordinates, the components of the velocity, and the time. 
However, in problems concerning the motion of electrons, and for that 
matter in many other physical problems also, we are usually concerned with 
forces of a somewhat special kind. Throughout the remainder of the article 
we shall assume that the force belongs to this special class. 

We consider four given functions of the coordinates and time, namely 

Vixi, Xo, Xs, t), ^l„(.ri, a-a, .Tg, /), (« ^ 1, 2, 3), 
and we assume that the components of the force are given by the formulae 



dxi dt ' [_ dxi dx^ . 

dX2 dt \_dx<i dXi 



^ ra£i _ dAi\ 



. :dA2 dAi~\ ,,, 



aF QA, . 



|_ 5^:3 dxi] ' |_ dx2 dXi J " 



Let us suppose, for purposes of illustration, that we are considering the 
motion of an electron. Then the physical interpretation of our assumption 

« L. A. MacColl, American Maihemalkal Monthly, Vol. 45 (1938), pp. 669-676. 
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concerning the force is the following. V{xi, X2, xs, t) is the potential energy 
of the electron in an electromagnetic field; that is 

V{xi, X2, Xz, t) = ~ &p(xi, X2, X3, t), 

where e is the absolute value of the electronic charge, and ip(xi, X2, X3, t) is 
the scalar potential of the field. The functions An(^'i, X2, a-'a, t) are related 
to the- components o„(a'i, X2, X3, t) of the vector potential of the field by the 
equations 

^^(.'^"1, X2, xs, t) = — ea„(xi, X2, Xs, t). 

The terms —dAn/dt are — c times the contributions of the vector potential 
to the components of the electric force. The quantity dAz/dx2 — dA^/dx^ 
is —tBi, where Bi is the a'l-component of the magnetic induction; and 
similarly for the quantities dA\/dxs — bAzjdxx and bA^lbxx — dAi/dx2.* 
In other cases also, equations (4), which may degenerate considerably, can 
be inteipreted without difficulty. 

Now we define a function L(xi, .tj, X3, Xi, ±2, X3, t) of the coordinates, the 
components of the velocity, and the time, as follows: 

L= - moc2(l - ^^c-2)'" - K + XiAi + 0-2^2 + XsA^, (5) 

We call this the Lagrangian function. 
We write the equations 

i~~l^-0. („ = 1,2,3), (6) 

at dXn ax„ 

carry out the indicated differentiations, and readily verify that the resulting 
equations are identical with those obtained by substituting the expressions 
(4) in equations (1). Hence, equations (6) are merely a form of the differ- 
ential equations of motion. We call equations (6) the Lagrangian equations. 
The chief importance of these equations is due to the ease with which they 
enable us to use coordinate systems which are not rectangular. This will 
be discussed in the final section. 

In the Newtonian case, i.e. the case in which the speed of the particle is 
small compared with the speed of light, the Lagrangian function reduces 
approximately to the form 

L= -moc' + ^ (x,' + xi 4- i/) - F + ±1 .4i + X2A2 + Xs As . (5') 

* These relations between the ^'s and the components of the vector potential, and 
between the partial derivatives of the A's and the components of the magnetic induction, 
are based upon the use of the M.K.S. system of units. If we measure the electromagnetic 
quantities in other units, certain constant proportionality factors may appear in the 
relations. 
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If we employ the function (5') in equations (6), we do indeed get the New- 
tonian differential equations. Since the constant term —ntoc^ is of no 
effect in the formation of the differential equations of motion, it is ordinarily 
omitted in writing the Newtonian form of the Lagrangian function. 

IV. HAAnLTON's Canonical Equations 
Let us write 

Pn-\- A„ = TTn. (7) 

Solving equations (2) for Xi, ±2, ±3, we get the result 

.T„ ^ cpAmo'c^ + pi" + p.' + ^V" 

(8) 
= c{wn - An)[mic' + (tti - Aif + {tt. - A,f + (tts - Az)V". 

Also, it is readily seen that the differential equations (1) can be written, 
with the aid of equations (7) and (8), in the form 



ox„ ox„ dx„ dXn 

= - ^ - c ^ KV + (x. - /Ii)^ + (x, - A,Y + (7r3 - A.fr. 



(9) 



Now let us define a function H(.\\, Xn, .vg, xi, X2, X3, t) as follows: 

H = dwoV + (xi - A,)- + (jra - A^f + (xg - .43)']'" + V. (10) 
Then equations (8) take the forms 

dH 



X. = ^ , (11) 



and equations (9) take the forms 

dll 

The function // is called the Hamiltonian function. The six equations 

(11) and (12), which are equivalent to the three equations (1), are called 
Hamilton's canonical equations of motion. These equations are of great 
importance in all of the deeper theoretical work in dynamics. 
An easy calculation shows that we have the identity 

II -\- L = TTiii + X2X2 + TVzXi. (13) 
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In the Newtonian case the Hamiltonian function given by (10) reduces 
approximately to the form 

H = moc' + ,— [(tti - Arf + (tt^ - A^f + (tts - A.f] + V. (10') 
Into 

The equations (11) and (12), with H given by (10'), are equivalent to the 
Newtonian differential equations of motion (!')■ Here again the constant 
term moc- is of no effect, and it is ordinarily omitted in writing the New- 
tonian form of the function H. The Newtonian forms of the functions H 
and L satisfy the identity (13), whether or not the constant terms moc^ and 
— moc^ are included. 

V. Static Fields of Force: The Energy Integral; Natural Families 

OF Trajectories 

By equations (11) and (12), we have the relation 

dt ^i Idx^ a7r„ ax„ dXni dt ' 

In particular, if no one of the functions V, Ai, A-i, A3 involves the time 
explicity, we have dH/dt = 0, so that the value of H remains constant 
during the motion of the particle. That is, under the condition stated we 
have 

moc'il - vc'T'" + Vixi, X2, X,) = constant. (IS) 

In the Newtonian case equation (15) reduces approximately to the form 

mc^ +~v -\- V(xi, XI, xz) = constant, 

which is equivalent to the equation 

^° v' + V(x,, X2, xs) = constant. (15') 

It is well known that this equation is a consequence of the Newtonian 
differential equations of motion. 

The left-hand member of equation (15') is the energy of the particle in 
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Newtonian dynamics, the first and second terms being the kinetic energy 
and the potential energy, respectively. The equation itself is called the 
energy integral/ Similarly, we call (15) the energy integral in relativistic 
dynamics, and we call the expression 

moc'[l-v' €-']-""+ V 

the relativistic energy. This energy is the sum of three parts: the proper 
energy mac^, the relativistic kinetic energy 

2r^ 2 — 2i~l/2 2 

woc [1 — V c ] — moc , 

and the potential energy V. 

The totality of possible trajectories of a particle in a static field of force 
forms a five-parameter family. We now see that if the field of force is 
static and of the kind we are considering now, the five-parameter family of 
curves consists of ^ i four-parameter subfamilies, each of which corresponds 
to a different value of the energy of the particle. Each of these four-param- 
eter subfamilies is called a natural family of trajectories. We proceed 
to derive the differential equations defining a natural family. 

If the constant in the right-hand member of equation (15) is denoted by 
the symbol E, we have the relation 

x,[l -H X,' -f xl'r - c[l - mlc\E - W'f^ (16) 

where 

X2 = dxs/dxi, xg = dxs/dxi. 
Hence, 

dt = c-'[l -\- x': + .v5Y'^[l - mlc* (E - VrV" dx,. 
From this, and the two equations 



d moX2 _ ^^ 1 • 

dt (1 - vH-'y- ~ ~d^2 ^' 



W0X3 dV , . 



\dx, dx,\ '''Idx, a^J' 

'L^ara dxi] ^'^\_dx2 dxsX 



dt (1 - ^2^-2)1/2 Q^^ 

it follows that we have the following system of differential equations defin- 
ing the natural family of trajectories corresponding to the total energy E: 

' In the theory of differential equations, an equation relating the unknowns involved 
in a system of differential equations, their derivatives of orders less than the highest orders 
appearing in the system, the independent variable, and one or more arbitrary constants, 
is called an integral of the system of differential equations. 
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,^_^ n, n.^w. d ( ,\{E- Vf - mWY\ 



2 411/2 



= i-{{E-Vf-moc'\ 



/ 2 / 

+ C[l + .T2 + 3:3 

[1 + x[" + :r 



/2._i/2 d ( I \{E — yf — nll^ 

3 J j~ I ^3 ~rr~; — '"2"; — '^ 

a^i\ L 1 + Xa + .T3 



2 4-11/2 
— niac 



(17) 



2 4il/2 



The equations which correspond to (17) in the Newtonian case are most 
readily obtained by going back to the Newtonian differential equations of 
motion and employing the integral 

An easy calculation, which is entirely parallel to the foregoing, gives us the 
following system of equations: 

+ [2„,i + .^' + .;'),- (4^ -|^]-[^-|^; 

+ [2^o(l + x. +X3)] ^ (^[ax^-aTj-^l^^^Jj- 
On comparing the systems of equations (17) and (17'), we get the follow- 
ing useful theorem. 

// the constants E, E*, mo, mo*, k, and the Junctions (of xi, Xn, X3,) V, Ai, Ai, 
A 3, V*, A 1*, Ai*, A 3* are such that we have identically 

(£- Vf -mlc' = k\E* - F*), 



(17') 



dA^ _dA2 ^ ^_ 

dx2 dx3 c{2mo*y'' 



p.4* dA*~\ 
[_ 6X2 dxs J ' 
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dxi 


dAx 

dX2 


dAi 


a.43 


dxz 


dxi 



c(2/Ho*)' 

k 



cilm^yi"- 



_ "a^ _ 6.4: 1 



Mi _ M*' 

3.T3 3a-i 



///c natural family of trajectories of a relativistic particle" {of rest-mass mo) 
moving with relativistic total energy E in the field of force derived from the 
functions V, Ai, A2, A3 is identical with the natural family of trajectories of a 
Newtonian particle {of mass nin*) moving with Newtonian energy E* in the 
field offeree derived from the functions V*, Ai*, A*, A*. 
In particular, the conditions of the theorem are satisfied if 

k = c(2wo)"', E* = c-2(2mo)->(£- - mlc'), m,* = m,, 

7 = K* = 0, A,* = A,, A,* = A2, As* = ^3. 

Hence, we have the corollary: 

In the case of an electrified particle moving in any static magnetic field the 
natural family of trajectories corresponding to any value of the energy given by 
relativistic {Newtonian) dynamics is identical with the natural family of 
trajectories corresponding to a certain other value of the energy given by New- 
Ionian {relativistic) dynamics. 

The equation 

E* = c~''^(2mo)~\E^ — mlc*) 

establishes a one-to-one correspondence between the physically significant 
{E ^ m-nc- and E* ^ 0) values of the relativistic energ>' E and the New- 
tonian energ}' E*. From this fact and the preceding corollary we get the 
following further result: 

fn the case of an electrified particle moving in any static magnetic field the 
total five-parameter family of trajectories given by relativistic dynamics is 
identical with that given by Newtonian dynamics. 

Of course, these peculiar properties of motion of an electrified particle 
moving in a static magnetic field are explained physically by the fact that 
the magnetic forces do no work, so that the speed of the particle, and 
consequently also its mass, remain constant during the motion. 

VI. Some Formulae from the Calculus of Variations 

This section is devoted to the derivation of some formulae from the 
Calculus of Variations which will be needed in the further discussion of the 

^I.e. a particle obeying the laws of relativistic dynamics. 
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dynamics of a particle. All constants, variables, and functions considered 
here are understood to be real. 

Let F{i, X, y, s, p, q, r) be a function of tJie seven arguments indicated, 
which, together with all of its partial derivatives of the first three orders, is 
continuous in a region R defined as follows: 

ai < i < Oi, 

hi < X < bs, 

R: Ci < y < C2, 

di < z < (h, 

p, g, and r unrestricted, 

the o's, 6's, c's, and (i's, being constants. 

Let x{t), y(t), z(t), <p{t), \l/{t), and co(/) be continuous functions with con- 
tinuous first derivatives, and let «, v, and 6 be parameters, independent of /, 
such that we have the relations 

61 < x{t) + e<p{t) < 62, 
ci < y{t) + #(0 < ^^2, {ai <t < (h). 

di < z{t) + 6oi{t) < d2, ... 

Let Ti and T2 be constants, and let h and ^2 be parameters, such that 
a, < Ti -{■ h < T2 -\- t2 < (i2. 
We now consider the integral 
/(€, V, 0, h, h) 



J T 



F{t, x-\- eip,y-\-r,^,z-\- do:, x' + €<p' , y' + t/i^', 2' + B^') dt. 

It can be shown without difficulty that the integral exists and is a differen- 
tiable function of e, rj, 8, h, h. We are interested in formulae giving the 
values of difdt, dl/dv, dl/de, dl/dh, dl/Bk at the point i= 7} = 6 = h = 
k = 0. 

' Since this section is purely mathematical, the constants, variables, and functions do 
not necessarily have any special physical significance. 

1" We treat the case of a function of seven arguments in order to fix the ideas, and 
because this is a case we shall meet in Section VII. However, the discussion applies 
essentially to other cases as well. In particular, in Section VII we shall also deal with a 
case in which F has only five arguments, z and r being absent. 
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By a well known theorem concerning the differentiation of definite 
integrals with respect to parameters,^^ we have 



dl 



Jri+t, L d{x-\-€ifi) d{x +€.^)J 

^ = -F[Ti-\- ti,x(Ti + h) -^ MT,-\- h), ■ ■ ■ ,^'iTi-\- k) -\- e<M'(Tx-\- Ix)], 
all 

^ = F[T2-\-h, x{T, + k) + ep(T2 + /2), ■ ■ • , z'{To + h) + eoi'{T2 + k)]. 
ah 

The formulae for 3//6ij and ai/dd are similar to that for dl/dt, and need 
not be written down. 

In particular, if [5//5e]o, etc. denote the values of the derivatives at the 
point e = Tj = 5 = /i = /2 = 0, we have 

[^1= -F\-Tux{T,), ■■■ ,2'(r01, 

The first three of equations (18) can be transformed to advantage, as 
follows. Integrating by parts, we obtain the formula 

and similar formulae for the integrals 

r\' ■^,F{i,x,---,z')dt 

Jti ay 



and 



j\'lF(,,..--,^)<U. 



'■ The theorem is given, often in the form of two separate theorems, in most works on 
Advanced Calculus and the Theory of Functions of Real Variables. See the bibliography. 
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It follows, therefore, that we have 

An important special case is that in which li and k are zero (so that the 
limits of integration are fixed), and 

^(^i) = <piT2} = rl^iTO = ^P{T2) = co(r:) = i.{T2) = 0. 
In this case we have in general 

/(., V, e, 0, 0) - 7(0, 0, 0, 0, 0) = e £" ,(0 g - I '^] dt 

r' ,.A^P d dFl _,, , , C' ,,JdF ddF~\,, 

+ o(e, V, 6), 

where o{€, tj, &) denotes a term, the exact form of which is unimportant, 
which is such that the expression 



approaches the limit zero as e, t], and 9 tend simultaneously toward zero. 

In particular, if the functions .v(/), y{t), s(/) satisfy the system of differen- 
tial equations ' 

ddF_dF^ ddF_3F^ ddF_dF^^ 

dtdx' dx ' dtdy' dy ' dl dz' dz ' ^ •' 

we have (for alt choices of the functions ip, ^, w subject to the conditions 
stated) 

/(e, V, e, 0, 0) - 7(0, 0, 0, 0, 0) - o{€, 71, d). (21) 
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Also, it can be shown without difficulty that in order that we have (21), for 
all such choices of ip, i/-, to, it is necessary that x, y, and z, satisfy the equations 

(20).* 
The last result can be stated in the following summary, and not quite 

explicit, form: If, and only if, the functions x(t), y(t), z(/) satisfy equations 
(20), the integral 

' F{t,x,y,z,x',y',z')dt (22) 



I 



is stationary with respect to infinitesimal variations of the functions x{t), 
y{i), z{t) which leave the terminal values unaltered. 

The problem of finding functions which render the values of definite 
integrals stationary is the chief subject of the Calculus of Variations. 

The equations (20) are called the Eulerian equations of the Calculus of 
Variations problem of making the value of the mtegral (22) stationary, or, 
as we usually say, of maximizing or minimizing the integral. 

VII. Hamilton's Principle and the Principle of Least Action 

We immediately recognize equations (6) as the Eulerian equations of a 
problem in the Calculus of Variations. Thus we have the following principle 
{Hamilton's principle): 

The particle moves, under forces of the type (4), so that the valtie of the integral 

Ldt, 



f 

Jt, 



with /i and In held fixed, is slationary with respect to infinitesimal variations of 
the functions x^it) which leave the initial and final points unaltered. 

The precise meaning of this is determined by the discussion given in 
Section VI. 

Hamilton's principle leads to the relativistic or Newtonian differential 
equations of motion, according as we use in it the function L given by (5) 
or by (5'). 

A little inspection suffices to show that the system of equations (17) is 
also the system of Eulerian equations of a problem in the Calculus of 

d dF dp 
* In brief, suppose that — — — — were not zero for some value of (. Then if we should 
^' dt dx' dx 

choose a function <p{l) which was (say) positive in the neighborhood of that value, and zero 
elsewhere, the integral 



!:>-{: ^m 



would have a value other than zero. We shall not give the actual proof here; it is to be 
found in the works on the Calculus of Variations cited in the bibhography. 
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Variations. Thus we get the so-called principle of least action, which can be 
stated as follows : 

The particle moves, in a static field of force of the type (4) , and with the 
prescribed total energy E, in such a curve that the value of the integral 

C'^ ([1 + x'^ + x:r[{E - v)\-' - mic'r + -4. + -4,:.; + ^3:^^) dx,, 

with the limits of integration held fixed, is stationary with respect to infinitesimal , 
variations of the trajectory which leave the end points unaltered. 

We have a precisely similar principle in Newtonian dynamics, but here 
the integral in question is 

f '"'' ([1 + x7 + x^'f"[2m,{E - V)r + ^1 + A2X2 4- A^x'z) dxu 

The last two integrals can be written more symmetrically, but not quite 
so explicitly, as follows: 

where Pi and ^"2 denote the end points of the trajectory, and ds^ = dxi + 
dxi -\- dxz. 

VIII. The Hamilton-Jacobt Theory 

Let us write 

W= r L[xi{t), X2{t), xsil), x[{t), X2{t), x,{t), t] dt. (23) 

We have already studied the variation of W when ti and t^ are held fixed, 
and the functions Xn{t) are varied in such a way that the terminal values are 
unaltered; and we have shown that under these circumstances the variation 
of W vanishes, to the first order of small quantities, in the natural motion. ^^ 
In the following we shall study the variation of W under some other 
conditions. 

Specifically, we shall study the quantity ATI' defined by equation (23) 
and the equation 

W + ^W = r ' L{xiit) + ?i(0, ■ ■ ■ , x[{t) + ^^{t), t] dt, 

" I.e. a motion satisfying equations (1). 
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where the functions x^{l) represent a natural motion, the ^„(/) are small 
functions, and Ah and A/2 are small parameters. 

It follows from the results of Section VI that we have (to within terms 
of the second order in small quantities)^^ 

AW = AkLMk), ■•• , X3(t2), h] 

■ - AhUxiih), •■■ ,X3{l{),h] 

t=-i tdx'nji^t^ 7.=! Lo3:„J(^ii 

= ■ AkUxiik), • ■ ■ , X3(t2), k] 
- AtiL[xi(h), ■■■ ,x'zik)>k] 

3 

+ E Wn(k)Uh) - Tr„{k)Ut>)]- 

n=l 

Let us write 

(Aa:„)ii = x,,{k + Ak) + Uh + A/.) - x^ik) - Uh) + x„{k) Ak, 

(Ax.)i = Xnih + Ak) + Uh + M) - x^ik) = Uh) + x'nOOA/i, 

so that (Aa;i)2, (Ax2)2, (Aa;s)2 are the coordinate differences of the termmal 
points of the varied and unvaried curves, and sunilarly (A%)i, (Aa:2)i, (Aa;3)i 
are the coordmate differences of the initial points. Then we have the 
formula 

AW = (Uxiik), ■■■] - Z) Mh)xUh)) Ak 

- (Uxyik), ■■■]-!: T„(k)xM) A/l 

3 

+ Z [ir„(k)i^x„)2 - 7rn(/i)(Aa-„)i], 

which, by equation (13), can be written in the form 
AW = -H[xi(k), ■■■]Ah-\- H[xi(k), ■ ■ ■] A/i 

+ E W.ikXAx^y. - Tn{k){AxM (24) 

Now, the integration ui (23) being taken over a natural motion of the 
particle, the value of W depends upon the initial instant, the initial coordj- 

" This is also the sense in which the following equations are to be understood. 
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nates, the initial components of velocity, and the final instant. It is 
necessary now to consider TT' as depending upon the following equivalent set 
of eight variables: the initial and final instants h and ^2, the coordinates 
(^11, ^21, -Tsi) of the initial point, and the coordinates (.ris, .T22, .T32) of the final 
point. Regarding W in this manner, we at once obtain the following 
relations from equation (24) 

dW „ dW ,._, 

^r = -^-h ^— = Tt.2, (25) 

df2 dx„2 

dt-i dx„i 

where H2 denotes H[xi(l2), ■ • ■ , Tri(j2), ■ ■ ■ , ^2] and Hi denotes H[xi{li), ■ ■ ■ , 

Let us now consider the partial differential equation 

^ + H{xi ,X2,X3, dW/dxi dW/dx2 , dW/dx3 , i) = 0. (27) 

at 

The preceding work shows that the function W we have been considering 
(with .Til, .T21, .V31, ti regarded as parameters, and with the symbols .V12, .T22, 
3:32, ^2 replaced by xi, .T2, :V3, t respectively) is a particular solution of this 
equation. We shall show that the cotnplete solution of this equation 
possesses remarkable properties in connection with dynamical problems. 

The complete solution of equation (27) is a function of .Ti, X2, ft'a, t, and 
of four arbitrary constants, of which one is merely additive, and can be 
neglected for our purposes. Let the solution be written 

W = W{xi, X2, X3, t, ai, a2, ai), 

where the a's are the three essential arbitrary constants. 

We write the equations 

1^ = ^. . (28) 

where the /3's are further arbitrary constants. These equations implicitly 
determine the .r's as functions of t and the six arbitrary constants ai, ■ ■ ■ , /3s. 
We also write the equations 

T— = Xn . (29) 

These equations determine three functions Tn. of the .r's, the a's, and t. 
In virtue of equations (28), the ir's are ultimately functions of /, the a's, and 
the i3's. 
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There is no reason to foresee a priori that the functions xi{t, ai, ■ ■ ■ , /Sg), 
• • • , Tsit, oi, ■ ■ ■ , 03) determined in this way, by means of the complete 
solution of equation (27), satisfy the differential equations of motion (11) 
and (12). Nevertheless, they actually do satisfy those equations, as we 
proceed to show. 

By equations (28), we have the relations 



dt dan dt m=i dctn dx^ 



(30) 



On the other hand, by (27) and (29), we have'* 

= -— — - + H{Xl, -Ts, A-3, Tl, X", ITS, t) 

dan \_ ot J 



dan dt ,„=i dir^ da„ 



The determinant 



aV 


dai aa:i 


aV 


das a.ti 



a IK Y dH dW 
dandl „i=i aTTn, dandx„ 

d'W 



(31) 



dai dXi 



d'W 



(32) 



daz dxs 
is not zero. For if it were, we would have a relation of the form 

^[air/a.vi, aTiya.va, aTTya.vs, .ri, :V2, :V3, t] = o, 

independent of the as. Now equation (32) is obviously distinct from (27), 
since it does not involve dW/dt. Hence, the vanishing of the determinant 
would imply that the function TJ'Cvi, x^, .rg, /, ai, a2, as) satisfies two distinct 
partial differential equations of the first order. This, however, is impossible 
when ir is the complete solution of (27); for an essential part of the concept 
of the complete solution of a differential equation is that the elimination of 
the arbitrary constants, from the solution and the equations obtained by 
differentiation, shall result in the given differential equation and no other. 
It follows, therefore, from (30) and (31) that 

dH 



Xtn ^~ 



aT„ 



We also have, by (29), 

aV 



Tn = 



dx„di 



+ E 



dW 



i dx^ aa;„ 



x„ = 



aV 



+ E 



aV dH 



dx„ dt „i=i dXm dXn air. 



{32,) 



'■' Since the function TI''(.Vi, x^., .ra, t, aj, aj, a^) satisfies equation (27) identically in the 
x's, t, and the a's. This remark apphes also in the case of equation (34), 
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On the other hand, we have 

■ 0-~r"+ffl = — + — + i:— ^'^ (34) 

By (33) and (34), we have the second set .of canonical equations 

■ _ ^^ 

TTn — r — - 

dx„ 

This completes the demonstration. 
If H does not involve the time explicitly, we can write 

W = S - El, (35) 

where E is an arbitrary parameter, and 5 is a solution of the differential 
equation 

H[xi, X2, X3, dS/dxi, dS/dx2, dS/dx^] = E. (36) 

The complete solution of (36) contains three arbitrary constants (besides 
the parameter E), of which one is merely additive, and can be neglected. 
It is easily seen that the solution of the canonical equations determined in 
the way described above, using the function W given by (35), and treating E 
as one of the a's, represents a motion of the particle with the total energy E. 
All of this theory holds both for the relativistic case and for the New- 
tonian case, the only difference being in the forms of the differential 
equations (27) and (36) in the two cases. 

IX. Curvilinear Coordinates 

In all of the foregoing we have employed rectangular coordinates, because 
they afford the simplest and most direct expression of the basic physical 
facts. However, in the solution of particular problems it is often more 
convenient to use other systems of coordinates. For this reason, we shall 
now formulate the more important equations in terms of general curvilinear 
coordinates. In this work, as in all work with general coordinate systems, 
we shall encounter concepts and relations which can be handled most 
perspicuously by means of the modem tensor calculus. Actually, the 
amount of tensor calculus we shall use is very slight, and no extended pre- 
limmary discussion is necessary in order to make the formulae intelligible. 
It will suffice to give occasional explanations of the notation, and of some 
of the concepts, as we proceed. Further information is to be found in the 
works cited in the bibliography. 

First consider the Lagrangian equations, which, as we have seen, are 
merely the Eulerian equations which follow from Hamilton's principle. 
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Now Hamilton's principle expresses a fact concerning the motion of a 
particle which is, by its very nature, independent of the choice of coordinates. 
Hence, the Lagrangian equations (6) hold in any coordinate system. How- 
ever, the form of the function L depends upon the particular coordinate 
sj'stem, and we must discuss the change of the form of the function resulting 
from a transformation of the coordinate system. 

In accordance with the common practice in the tensor calculus, we shall 
now denote the coordinates by the symbols x , x , x- , instead of by the 
symbols Xi, X2, X3. 

In rectangular coordinates the differential distance ds between the points 
{x^, x^, x^) and (x^ + dx^, x"^ + dx"^, x^ + dx^) is given by the shnple formula 

ds^ = dx^^ + dx^' + dx^^ , 

but this is highly special; in general coordinates we have 

ds' = t tgUx,x\x')dx-dx\ (37) 

where the g's are functions which depend upon the particular coordinate 
system under consideration. It is understood that gmn = gnm ■ Hence- 
forth, we shall write (37) in the form 

ds^ = g„„(.-vS x', x^)dx^dx^, (38) 

and we shall observe this general rule throughout: When the same literal 
index occurs twice in a term, once as a subscript and once as a superscript, 
that term is understood to be summed for the three values of the index. 
We now have the result 

v^ = [ds/dty- = g„.(x\ x\ x^)x"-z^, 

and 

nw'il - v-^tr-^y^ = moc'{l - (r-''g„„x'-x-Y'\ 

The function V{x^, x-, x^, t) is a scalar. That is to say, when the coordi- 
nate system is changed, the first three arguments of the function are replaced 
by their expressions in terms of the new coordinates, and so we obtain a 
function which is of a new analytical form, but which has the same value 
as the original fimction at each point of space. 

Now we consider the expression 

Aix^ -\- A^x^ -\- A3X\ 

In rectangular coordinates this is the scalar product of the vectors (^1, 
Ai, A^ and (x', x^, x^). The expression retams its form and interpretation 
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under changes of the coordinate system, provided (as the notation implies) 
{A-i, Az, As) is treated as a covariant vector.^^ 

With these understandings as to the significance of the symbolism, we 
can now write down the following general expressions for the Lagrangian 
function L in the relativistic and Newtonian cases, respectively, 

T- 2 r< —2 -m .nil/2 rr i i -tn 

L = -mac \\ — c gmnX X \' - V + A^x , 

= —moC + J gmnX .X — V -\- A^^x . 

These hold for any coordinate system; and from the appropriate one of these, 
and the Lagrangian equations 

ddL^dL 

dt dx" ar" ' 

we obtain the relativistic or Newtonian differential equations of motion in 
any coordinates. 

Now let us consider the Hamiltonian canonical equations. 

We have already agreed to consider (^i, A^, A^) as a covariant vector. 
We now make the same convention in regard to (tti, iri, tts). Then it read- 
ily follows that the equations 

a^» - - («) 

^^ Suppose that with a point P (which may be either a special point or a typical point), 
and with each coordinate system, we have associated an ordered triple of numbers. 

If the triples of number (ffli, Oa, Oa) and (ffli', a^', oa') associated, respectively, with any 
two coordinate systems {x^, x^, x^) and x^', x'^' , x^') satisfy the relations 

dm' = r— o«, 
o.t" 

the numbers (tii, 02, os) are said to be the components of a covariant vector in the coordinate 
system {x^, -i*', x^). (It is understood, of course, that the partial derivatives are evaluated 
at the point P.) 

On the other hand, if the triples of numbers {a', a^, a?) and (a^', a^', a^') associated with 
the typical coordinate systems (.t', x^, x^) and («'', x'^' , x^') satisfy the relations 





dx'^' 


a"* 


= a", 




dx" ' 



the numbers (a^, 1^, 1^) are said to be the components of a contravariant vector in the 
coordinate system (x^, x'^, x'). 

These concepts agree only in part with the ones used in the elementary theory of 
vectors. From our present standpoint, the only vectors used in the elementary theories 
are those which are defined with reference to rectangular coordinate systems. When other 
coordinate systems are used (e.g. cyhndrical coordinates), the vectors, defined in terms of 
rectangular coordinates, are merely resolved along the tangents to the coordinate curves. 
The components obtained in this way are not the same as the components considered in 
the tensor calculus, which we are using here. 
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are tensor equations; and since they hold when the coordinates are rec- 
tangular, they hold for all coordinate systems/^ 

We let g™" denote g~^ times the cofactor of the element g^n in the deter- 
minant 



gn 


gll 


gU 


^21 


gl2 


g23 


g31 


g32 


^33 



Now we write 

H = c[mlc' + ^"{Tr. - ^„,)(7r„ - A.)Y" + V (41) 

for the relativistic case, and 

H = moc^ + (2mo)-'r"(vr.. - ^l™)(7r„ - A„) + V (41') 

for the Newtonian case. We see that these expressions specialize into the 
ones given earlier for the Hamiltonian function when the coordinates are 
rectangular. 
E, L, and Xni" are all scalars. Consequently, the equation 

H -\- L= TT^i" (42) 

is a tensor relation; and since it holds when the coordinates are rectangular, 
it holds for all coordinate systems. 
The Lagrangian equations can be written in the form 

Let us consider the variation of the function L resulting from small 
variations of the x's and x's. By (40) and (43), we have the relation 

arc" a.T" 

= 5(7r„i") + (t„5a:" - x"57r„). 

It follows from (42) and (44) that the variation of H resulting from small 
variations of the x'& and the tt's is given by the formula 

SH = x"bir„ — TTr^dx". 

^* The argument is explained in detail in the works on the tensor calculus cited in the 
bibliography. 
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From this it follows that we have the Hamiltonian canonical equations 
.„ _dll . _ _dH 

dTn ™" 

in any coordinate system. 

We have already seen how to state Hamilton's principle in terms of general 
coordinates. 

In the relativistic case the principle of least action takes the form: The 
particle moves, in a static Ji eld of force of the type (4), and with the prescribed 
total energy E, in such a curve thai the value of the integral 



/ 



/r dx" dx"~y^\.^ -,.s2 -2 2 2,1/2 , , ^^M ^ 1 



with the limits of integration held fixed, is stationary with respect to infinitesimal 
variations of the trajectory which leave the end points unaltered. The corre- 
sponding form of the principle for the Newtonian case is obvious. 

We are now in a position to dispose very quickly of the problem of formu- 
lating the Hamilton-Jacobi theory in terms of general curvilinear 
coordinates. 

The general form of the Hamiltonian function being given by (41) (for 
the relativistic case) or (41') (for the Newtonian case), we can at once write 
down the partial differential equation 



dw 

Let 



— + I-I{x\ x\ x\ dW/dx\ dW/dx\ dW/dx\ t) = 0. (45) 

dt 



W = W{x\ x\ x', t, a\ a^ a^) 



represent the complete solution of (45), without the irrelevant additive 
constant of integration. 

Our chief problem is that of proving that the functions x"{t, a\ a^, «■', 
ft , /32 , 0s), TTnit, «S a^ a^ ft , /?2 , &z) determined by the equations 

dW ^ dW 

a^ - '^" ' ^ - ''" ' 

where the /3's are further arbitrary constants, satisfy the canonical equations 

.„ dll . _ dH 

^ — 1. — J Tr„ — — — - . 

dir„ ox" 

Now, referring to the proof given in Section VIII for the special case of 
rectangular coordinates, we see at once that nothing in the proof depends 
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upon the special forms which the Hamiltonian function and equation (45) 
assume in those coordinates. Hence the proof already given applies imme- 
diately to the present general case. 

Similar remarks apply also to the case in which H does not involve the 
lime explicitly, and in which we write 

ir = 5 - Et, 

where S is the complete solution (without the additive arbitrary constant) 
of the equation 

H{x\ .r^ x\ dS/dx\ dS/dx\ dS/dx') = E. 
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